Coulomb deexcitation of pionic hydrogen within close-coupling method 
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The Coulomb deexcitation of (np) atom in collisions with the hydrogen atom has been studied 
in the fully quantum mechanical close-coupling approach for the first time. The total Coulomb 
deexcitation cross-sections of nl — > n'l' transitions and Z-averaged cross-sections are calculated for 
n = 3 — 8 and at the relative energies E = 0.01 — 100 eV. The strong interaction and vacuum 
polarization shifts of ns-states are taken into account. It is shown that the An > 1 transitions are 
very important and make up a substantial fraction of the Coulomb deexcitation cross-section (up 
to ~ 47%). 
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Exotic hydrogen atoms are formed in highly excited 
states and the evolution of their distributions on the 
quantum numbers and kinetic energy is defined by the 
competition of both coUisional and radiative deexcita- 
tion processes during the atomic cascade. The interest in 
the problem attracted special attention after the exper- 
imental observation of a high energy fraction (3> 1 eV) 
in the energy distribution of the pionic atoms at the in- 
stant of the charge exchange reaction ir^p ir^n was 
obtained P,|3j0;IJI- The observed Doppler broadening of 
the neutron time-of-flight (NTOF) spectra in the exper- 
iments are attributed to the Coulomb deexcitation (CD) 
process, where the energy of the transition is shared be- 
tween the colliding objects. In order to analyze the above 
mentioned and the new experiments 0, S S 13 on the 
precision spectroscopy of K X-rays from pionic atoms, 
the kinetic-energy distribution of pionic hydrogen must 
be calculated from a realistic cascade model. 

The first theoretical study has been carried out in the 
framework of the two-state semiclassical model with 
some additional approximations (e.g., n 3). Later, 
the process was considered within the advanced adia- 
batic ^3 (see also references therein), and the classical- 
trajectory Monte Carlo ^3 approaches. While for the 
low n states the approaches [3, ^3 cannot be expected 
to give a reliable description of the CD process, the CD 
cross sections obtained in for the n—i-n—1 transitions 
(n = 3 — 5) are smaller about one order of the magnitude 
than in |9( and too small to explain the experimental 
data 1^3- Therefore, in the most important region of 
n = 3 — 7 relevant for the atomic cascade the theoretical 
results are rather undefined and this process is regarded 
as the least known process of the pionic cascade. 

In this paper, we report the results for the CD process 
of 7r~p atom in collisions with H obtained for the first 
time in the framework of the fully quantum mechanical 
approach. We use the close coupling (CC) method, in 
which the scattering processes 



(tT p)nl + His (tT p)n'l' + His 



(1) 



such as elastic scattering {n' = n, I' — I), Stark transitions 
(n' = n, I' y^l), and CD (n' < n) are described in unified 
manner. The method has been recently employed for the 
study of the CD process in muonic hydrogen [Tjl , where 
substantially new results in comparison with the known 
in the literature |^, ^3 (see references therein) have been 
obtained. Unless otherwise stated the atomic units {h = 
e — niemp/ (me + mp) = 1) will be used throughout the 
paper. 



The Hamiltonian of the system ((tt p)ni + His) (after 
separation of the cm. motion) is given by 



i/ = - — Ar 
2 m 



-/ie.(p) + /^ff(r) + F(r,p,R), (2) 

where m is the reduced mass of the system, R is the ra- 
dius vector between the cm. of the colliding atoms, p 
and r are their inner coordinates. The interaction po- 
tential, l/(r,p, R), is a sum of four Coulomb pair inter- 
actions between the projectile atom and the target atom 
particles, h^xip) and /i//(r) are the hydrogen-like Hamil- 
tonians of the free exotic and hydrogen atom, whose 
eigenfunctions together with the angular wave function 
YlaCR-) of the relative motion form the basis states, 
\ls,nl,L : JM), with the conserving total angular mo- 
mentum (JM) and parity tt = (—1)'+^. The total wave 
function of the system can be expanded as follows 

^■'/'^r, p,R)^R-'J2 GilLim^s, nl, L : JM). (3) 

nlL 

The expansion (3) leads to the close-coupling second 
order differential equations for the radial functions of the 
relative motion, G;{^(i?), 
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{R)GiTi,AR)- (4) 



The channel wave number is defined as k"^ = 2m{Ecm + 
Enolo - Eni), where Ecm and [uqIq) are the energy of the 
relative motion and the exotic atom quantum numbers 
in the entrance channel, respectively. The bound energy 
of 7rp atom, Eni = Sni + Ae„i , includes the eigenvalueof 
hex{p), £nh and the energy shift, Ae„; = Ae^^ , 
due to the vacuum polarization and strong interaction. 
Hereafter, the energy Ecm will be referred to e„;/o in the 
entrance channel (we assume that Ae„/^o = 0)^. 
The matrix elements of the interaction potential. 



I' L' .nlL 



{R) = {Is, n'l', L' : J M\V\ls,nl, L : JM) 



(5) 



The strong interaction shifts are calculated acco rding to Ae^*J = 
Aef/7n3, where Aei^ = -7.11 eV from E3)- The next val- 
ues (-3.24 eV (Is), -0.37 eV (2s), and -0.11 eV (3s)) for 
the vacuum polarization shifts are used in the calculations. For 



n > 3 we assume A£„ 



Ae" 



2 



are obtained by averaging it over the electron wave func- 
tion of the Is-state and then applying the multipole ex- 
pansion. The integration over (p, R) reduces the matrix 
elements (5) to the multiple finite sum. 

At fixed Ecm the coupled differential equations (4) for 
the given J and n values are solved numerically by Nu- 
merov method with the standing-wave boundary condi- 
tions involving the real symmetrical ii'-matrix related to 
T-matrix by the equation T = 2K/{I — iK). All exotic 
atom states with n values from 1 up to hq have been in- 
cluded in the close-coupling calculations. The following 
cross sections will be discussed: 

- partial cross section 

- total cross section of the transition nl n'V 

cT,,i^,,i'{E)=Y^aii^^,i,{E), (7) 
J 

- total cross section from ns - state 

(Tns^r.'{E)=Y,(^ns^n'v{E), (8) 
I' 

and the Z-averaged cross section [l > 0) 

ai>?(i?) = Yl + l)ani^n'i'{E). (9) 

l>0,l' 

We assume the pionic atom states with / 7^ are degen- 
erated and so the /-averaged cross section (9) can be used 
to illustrate the basic features of the CD process. 

The CC calculations of the CD cross sections Eqs. (6- 
9) were carried out for {tt^p)„i+H collisions for n = 3 — 8 
and at energy range Ecm — (0.01 — 100) eV. Summation 
over the partial waves in Eq. (6) was done up to a value 
Jmax until an accuracy better than 0.1% was reached at 
all energies. The analysis of the J dependence of the par- 
tial cross sections (t;[;^„,;, shows that most the CD cross 
sections give the partial waves with rather low J values as 
compared with the elastic scattering and Stark transition 
processes. Besides, our study shows that a significant un- 
derestimation of the CD cross sections (up to one order of 
the magnitude) can be obtained if one neglects the higher 
multipole terms in the expansion of the coupling matrix 
(5). Such a strong effect together with the relatively low 
J values involved in the inelastic transitions are related 
with the fact that the CD process is determined by the 
comparably more short-range interaction. 

In order to illustrate the influence of the ns state en- 
ergy shifts on the CD cross sections, the calculations were 
performed both with and without taking energy shifts 
into account. The effect is the most pronounced for the 
low-lying states and some of the results for n = 3, 4 are 
shown in Figs. [5] and presented in Table I. 

The CD cross sections for a number of transitions 
3/ 21' with or without the energy shifts of the ns-states 
are shown in Fig. ^ In the case A£„s ^ the energy 
dependence of the cross sections changes crucially. The 
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FIG. 1: The CD cross sections (7^1^21' for (7rp)3; -\- H collisions 
calculated both with (thick lines) and without (thin lines) 
taking the ns-state energy shifts into account. 

cross sections corresponding to the 3p(3d) — > 2s, 2p tran- 
sitions at energies above and below the threshold of the 
Stark transitions 3s 3p, 3d are strongly suppressed. 
The effect is more clearly observed at energies in the 
vicinity of 0.1 eV (the suppression is more or about one 
order of the magnitude). At Ecm ^ |Ae„s| the differ- 
ences between the corresponding pair of the curves (see 
Fig. ^1 are negligible at Ecm > 15 cV for n = 3 and at 
Ecm > 1 eV for n = 4. On the contrary, the transitions 
3s — > 2s, 2p are enhanced at the threshold by a factor 
about 3 — 4. The main reason of such behavior of the 
CD cross section is related with the inelastic nature of 
the Stark transitions between the splitted upper states. 
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FIG. 2: The /-averaged CD cross sections a'^^ (thick lines) 
and a^^2 (thin lines) for {tvp)4 + H collisions calculated with 
(solid lines) and without (dashed lines) taking the s-states 
energy shifts into account. The CD cross sections cr4s_3 cal- 
culated both with (double-dot-dashed line) and without (dot- 
dashed line) taking the s-state energy shifts into account are 
also shown. 

The similar regularities are also observed in the l- 
averaged CD cross sections. In Fig. |21(see also Table 1) 
the /-averaged CD cross sections for the 4^3 and 4^2 
transitions and total cross sections from 4s-state are pre- 
sented. One can see that the maximal suppression due to 
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TABLE I: The CD cross sections crl^^, for {tvp)„ + H collisions calculated in quantum-mechanical CC approach. 
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a"32 

(Arf)) 
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0.409 
2.261 


0.231 
1.315 


0.289 
1.293 


0.145 
1.303 


0.138 
1.275 


0.143 
1.711 


0.139 
0.770 


0.102 
0.293 


0.072 
0.148 


0.046 
0.075 


0.022 
0.030 


0.017 
0.021 


0.012 
0.015 


0.007 
0.008 


0.003 
0.004 


0.002 
0.002 


0'43 

(A^) 

0"42 
(A^) 
"42 


1.342 
2.558 
0.438 
0.977 


1.076 
1.728 
0.491 
0.946 


0.858 
1.402 
0.489 
0.942 


0.638 
1.052 
0.404 
0.817 


0.958 
1.004 
0.387 
0.721 


0.642 
0.969 
0.252 
0.456 


0.599 
0.814 
0.161 
0.250 


0.399 
0.483 
0.100 
0.129 


0.245 
0.280 
0.060 
0.071 


0.132 
0.145 
0.032 
0.036 


0.058 
0.062 
0.014 
0.015 


0.042 
0.044 
0.010 
0.011 


0.030 
0.031 
0.007 
0.008 


0.017 
0.017 
0.004 
0.004 


0.009 
0.009 
0.002 
0.002 


0.006 
0.006 
0.001 
0.001 


0"54 
0"53 
0"52 


3.497 
0.929 
0.188 


2.672 
1.118 
0.292 


1.750 
0.677 
0.199 


1.201 
0.417 
0.122 


1.055 
0.355 
0.103 


0.769 
0.251 
0.074 


0.531 
0.177 
0.069 


0.314 
0.091 
0.040 


0.202 
0.050 
0.022 


0.125 
0.025 
0.011 


0.060 
0.010 
0.005 


0.044 
0.007 
0.004 


0.033 
0.005 
0.003 


0.020 
0.003 
0.001 


0.013 
0.002 
0.001 


0.012 
0.001 
0.001 


0"65 
0"64 
0"63 


3.644 
1.446 
0.436 


2.510 
1.309 
0.349 


2.423 
1.232 
0.265 


1.544 
0.903 
0.171 


1.121 
0.683 
0.115 


1.012 
0.603 
0.098 


0.636 
0.377 
0.065 


0.298 
0.217 
0.032 


0.169 
0.129 
0.017 


0.098 
0.068 
0.008 


0.041 
0.028 
0.003 


0.031 
0.021 
0.002 


0.024 
0.015 
0.002 


0.016 
0.009 
0.001 


0.015 
0.005 
0.001 


0.020 
0.004 
0.001 


0"76 
0"75 
0"74 


5.647 
1.545 
0.300 


4.321 
1.272 
0.280 


3.202 
0.793 
0.154 


2.636 
0.696 
0.156 


2.157 
0.762 
0.162 


1.747 
0.527 
0.114 


1.247 
0.319 
0.069 


0.590 
0.150 
0.037 


0.346 
0.095 
0.019 


0.190 
0.053 
0.010 


0.082 
0.023 
0.004 


0.061 
0.018 
0.003 


0.045 
0.014 
0.002 


0.027 
0.009 
0.002 


0.023 
0.007 
0.001 


0.038 
0.008 
0.002 


0"87 
0"86 
0"85 


7.480 
1.340 
0.571 


4.541 
0.925 
0.340 


4.095 
0.822 
0.315 


3.379 
0.719 
0.362 


3.217 
0.624 
0.234 


2.379 
0.485 
0.222 


1.757 
0.378 
0.170 


1.003 
0.190 
0.078 


0.662 
0.130 
0.044 


0.371 
0.072 
0.023 


0.159 
0.031 
0.010 


0.118 
0.023 
0.008 


0.086 
0.017 
0.006 


0.051 
0.011 
0.004 


0.042 
0.010 
0.003 


0.077 
0.016 
0.004 



the energy shift of 4s state is weaker and less or about two 
times at very low energy both for 4 — s- 3 and 4 — > 2 tran- 
sitions, while at Ecm > 1 eV does not exceed 15%. The 
important observation follows from Fig. El the CD cross 
sections for the Art — 2 transition are strongly enhanced, 
in contrast to the calculations 0, ^ , and their relative 
contribution to the total CD transitions 4 — > n' < 3 
varies from 20% up to 39%. 
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FIG. 3: The i-averaged CD cross sections (t'j''°_j^ for {■Kp)„+H 
collisions with n = 5 (solid line), 6 (dashed), 7 (dot-dashed) 
and 8 (dotted). The results of the parameterization |1H| for 
n — 5 and 7 (thin lines). 

Beginning from the paper Q , it is commonly believed 
that the CD cross sections for An = 1 transitions behave 
like rf /Ecm- The n and Ecm dependences of the l- 
averaged CD cross sections (/ > 0) for n = 5 — 8 are 
shown in Fig. |3 for the An=l transitions. The straight 
lines show the results of the parameterization 1 1] , based 



on the CD cross sections |^ for n < 7. As it is seen 
from Fig. O the energy dependence of the CD cross sec- 
tions only in the region 1 ;$ -E ^ 10 eV, as a whole, 
is in a qualitative agreement with the one like 1/Ecm- 
But at the energies beyond this interval the CC results 
don't confirm the 1/E energy dependence of the cross 
sections. In particular, at < 1 eV and n>4 our results 
show weaker and rather 1 / ^/E behavior of the CD cross 
sections. 

Concernig the power n-dependence nearly to n'^ with 
7 > 2, the present consideration doesn't confirm that the 
CD cross sections have such a scale factor depending on 
n (see Fig. O and Table 1). Moreover, for n = 4 6 the 
non-monotone behaviour of the as a function of n 
is revealed. 
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FIG. 4; The same as in Fig. 3 but for transitions n n — 2. 

The present calculations show that An = 1 transitions 
dominate. According to our CC calculations (see Fig. 0| 
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FIG. 5: The CD cross sections from s-states a„s^n-i for 
n = 5 — 8. The notations of curves are the same as in Figs. 3 
and 4. 

^ We thank Dr. T. Jensen for the preliminary results of the anal- 
ysis with the present CD cross sections of the NTOF spectra . 
Our results explain the high energy component around 105 eV 



and Table I) , the An > 1 transitions make up a substan- 
tial fraction of the total CD cross sections (20-47 %) for 
n > 4 at all the energies under consideration. It is worth- 
while noting, that there is no simple n dependence of the 
CD cross sections for the An > 1 transitions. The reg- 
ular n dependence appears only at the relative energies 
more or about 50-70 eV. 

In summary, the CD process in pionic hydrogen has 
been studied for the first time in the fully quantum me- 
chanical CC approach, involving simultaneously the elas- 
tic and Stark processes. The obtained results reveal the 
new knowledge about the CD process and are very im- 
portant for the reliable analysis of the K X-ray yields 
and NTOF spectra.^ 

We are grateful to Prof. G. Korenman for fruitful dis- 
cussions. This work was partially supported by Russian 
Foundation for Basic Research, grant No. 03-02-16616. 

due to the 5 — > 3 transition and lead to a very good agreement 
with the experimental weight of the 3^2 component at 209 eV. 
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